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Testing black hole superradiance with pulsar companions 
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We show that the magnetic dipole and gravitational radiation emitted by a pulsar can undergo 
superradiant scattering off a spinning black hole companion. We find that the relative amount of 
superradiant modes in the radiation depends on the pulsar’s angular position relative to the black 
hole’s equatorial plane. In particular, when the pulsar and black hole spins are aligned, superradiant 
modes are dominant at large angles, leading to an amplification of the pulsar’s luminosity, whereas 
for small angles the radiation is dominantly composed of non-superradiant modes and the signal is 
attenuated. This results in a characteristic orbital modulation of the pulsar’s luminosity, up to the 
percent level within our approximations, which may potentially yield a signature of superradiant 
scattering in astrophysical black holes and hence an important test of general relativity. 
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It has been known since the 1970s that low-frequency 
waves can be amplified upon scattering off a rotating 
black hole (BH) [lHa]. Superradiant scattering occurs for 
waves of frequency ui < mfl, where f1 denotes the angular 
velocity of the BH horizon and m the azimuthal number 
that characterizes the wave’s angular momentum. Such 
waves carry negative energy and spin into the BH from 
the point of view of a distant observer, which is possible 
inside the ergoregion that surrounds a rotating BH. Such 
an observer thus sees the BH lose a fraction of its mass 
and angular momentum in the scattering process (see @ 
for a recent review of this topic). 

Superradiance is thus inherently associated to the 
space-time distortion produced by rotating BHs. Finding 
observational evidence for this process would then pro¬ 
vide an extremely important test of Einstein’s general 
relativity and a unique insight into the properties of the 
most compact objects in the universe. 

Pulsar-BH binaries have long been regarded as the 
‘holy grail’ for testing the fundamental theory of grav¬ 
ity, combining one of the most accurate ‘clocks’ in the 
universe with the strong space-time warping produced 
by a BH. Finding these systems is, in fact, one of the 
key goals of the Square Kilometer Array [7|, as well as 
for current gravitational wave observatories such as Ad¬ 
vanced LIGO/VIRGO i, where detection of up to a few 
binary mergers per year may be expected [g] . 

In this Letter, we show that the radiation emitted by a 
pulsar can undergo superradiant scattering off a rotating 
BH companion, constituting the first example of an as¬ 
trophysical system where observational evidence for BH 
superradiance can be found. 

We begin by considering the angular velocity of a Kerr 
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BH of mass M and spin J = aMc, which is given by: 


n = 
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where r + = (GM/c 2 )( 1 + \/l — a 2 ) denotes the radial 
position of the outer horizon and a = Jc/GM 2 is the di¬ 
mensionless spin parameter. This implies that, for stellar 
mass BHs in the range 10 Mq — 30 Mq, superradiant am¬ 
plification can only be significant for waves below the 
1 kHz — 10 kHz frequency range, particularly since it is 
suppressed for large multipole numbers l and m 

Such low frequency radiation is naturally produced by 
pulsars. These rapidly spinning neutron stars harbor ex¬ 
tremely large magnetic fields formed in the collapse of 
the parent star. For ordinary pulsars, these can reach 
10 8 — 10 9 T, and even larger values for magnetars. This 
makes pulsars powerful emitters of magnetic dipole ra¬ 


diation at the rotation frequency ui p 


1 Hz - 10 kHz, 


the upper bound correspon ding to the rapidly spinning 
millisecond pulsars (see e.g. 10j). 

In addition, spinning neutron stars may also be an 
important source of gravitational waves due to devia¬ 
tions from axial symmetry. Several sources of asymmetry 
have been discussed in the literature, such as anisotropic 
stresses supported by the neutron star’s solid crust, a 
misalignment between the rotation axis and the princi¬ 
pal axis of inertia as a result of its violent formation, or 
a misalignment between the rotation and magnetic axes 
(see e.g. ;IH). Due to its quadrupole nature, gravita¬ 
tional radiation is generically emitted at both the pul¬ 
sar’s angular frequency and twice its value, in amounts 
depending on the structure of the deformation. 

It is thus clear that pulsars can emit both electro¬ 
magnetic (EM) and gravitational waves (GWs) with fre¬ 
quency below the threshold for superradiant scattering, 
potentially providing two different channels for astro- 
physical observations. The question that remains to be 
addressed is whether this radiation has the correct mul¬ 
tipolar character to extract energy and spin from a BH 
companion. 
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Let us first review the basic features of EM and GW 
scattering in the Kerr spacetime. A convenient frame¬ 
work to study these processes is the Newman-Penrose 
(NP) formalism [l2| , where one projects the Maxwell and 
Weyl tensors describing the waves along a null tetrad. 
For the Kerr spacetime, a useful choice is the Kinnersley 
tetrad [l3| , where two of the 4-vectors coincide with the 
principal null directions of the Weyl tensor. This allows 
one to construct complex scalar functions, T, s , known as 
the NP scalars and which encode the different spin com¬ 
ponents of the EM (s = ±1) and GW (s = ±2) radiation. 

These quantities satisfy independent scalar wave equa¬ 
tions that can be solved using separation of variables 0. 
In particular, they admit a generic mode expansion in 
Boyer-Lindquist coordinates (t, r, 9, (f>) 0 of the form: 

T s = J2 e~ iUlt e im *sSim(d)sRim(r) + (oj^- w)(2) 


For simplicity, we will first focus on modes of the form 
e~ lut and discuss the e lut modes below. The angu¬ 
lar functions s Si rn (9)e' lm<l> correspond to spin-weighted 
spheroidal harmonics, which reduce to the well-known 
spin-weighted spherical harmonics s Yi m (9,(j)) 0 for 
aco < 1. The problem then reduces to a single radial 
wave equation, known as the Teukolsky equation Q: 
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where, in units such that G = c = h = 1, K(r) = (r 2 + 
a 2 )ui — ma 1 A = r 2 — 2Mr+a 2 and the angular eigenvalue 
can be written as a series expansion s Ai m = 1(1 + 1) — 
s(s + 1) + Efc°l c fc(« w ) fc 0- 
Although no exact solution of this equation is known, 
approximate solutions can be found for low frequencies, 
u> -C r7 , in two overlapping regions. These “near” and 
“far” solutions, valid for r — r + -C w _1 , and r — r+ r + , 
respectively, can then be matched to produce a com¬ 
plete solution. We then obtain a relation between the 
incoming and outgoing radiation, and consequently the 
energy flowing into the BH horizon. This is encoded in 
the gain/loss factor Q: 

dEout/dt ^ dJ'Jii!! / di 

dEin J dt dEi n J dt 

+ . (4) 
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where vo = (to — mfi)r + , e = (r+ — r_)/2M and Ci s = 
(l + s)\(l — s)\/(21)1(21 + 1)!. We thus conclude that su¬ 
perradiant scattering, s Zi m > 0, occurs when oj < mil 
for co-rotating wave modes, m > 0. Also, amplification 
is larger for the lowest multipoles and, since l > |s|, it is 
more pronounced in the EM case (at low frequencies). 


We note that s Z; m (w) = ui) for arbitrary fre¬ 

quencies Q, so that superradiance will occur for m < 0 
for the wave modes of the form e lwt . 

Superradiant amplification is, however, much larger for 
higher frequencies, in particular for near-extremal BHs 
(a ~ 1). In this regime, we need to employ numeri¬ 
cal methods to solve the Teukolsky equation. A simple 
procedure, first used in Q|, consists in numerically inte¬ 
grating an ingoing wave solution at the horizon up to a 
large distance, where one can extract the coefficients of 
the incoming and outgoing radiation, as we describe in 
detail in a companion article m 

For GWs, we obtain a maximum amplification Z max ~ 
1.02 for the l = m = 2 mode with a = 0.999 and 
ui — 2fl. In the EM case, the maximum gain is consider¬ 
ably lower, Z max ~ 0.044, for l = m = 1, in agreement 
with [3j]. For non-superradiant modes, Z < 0 and \Z\ 
increases with the frequency. In particular, the lowest 
non-superradiant multipoles approach the maximal ab¬ 
sorption limit, Z ~ — 1, for u> < m£l. In addition, for a 
given frequency, higher multipoles exhibit progressively 
smaller gain/loss factors, such that only a finite num¬ 
ber of modes will effectively be relevant in a scattering 
problem (see also 0 ). 

A physical wave, described by real fields, generically 
contains both e ±lult frequency modes, as well as different 
(l,m) multipoles. Whether an overall amplification of 
the signal occurs upon scattering off a Kerr BH then de¬ 
pends on the relative amplitude of the different incoming 
modes, which we will now determine for the pulsar-BH 
binary. 

We consider the limit where the orbital distance is 
large, d>A>r + ,l? p , where A is the wavelength of the 
radiation and R p the neutron star’s radius. In this limit, 
it is a good approximation to consider the electromag¬ 
netic and gravitational fields generated by the rotating 
neutron star in flat space, which then give the incident 
wave for the scattering problem. 

Since in this limit the orbital period will largely exceed 
the pulsar’s rotational period, we may first study the 
scattering problem for a given pulsar position (d, 9 p: <fi p ) 
in the BH frame and then include the orbital variation 
of these parameters. 

It is also a sufficiently good approximation to treat 
the pulsar as a point-like spinning magnetic dipole and 
mass quadrupole. Here we will focus on binary systems 
where the pulsar and BH spins are either aligned or anti¬ 
aligned, although the same method can be applied to 
arbitrary configurations. 

For a misalignment angle a between the pulsar’s spin 
and its magnetic dipole moment, we have: 


1 „ 

M„ = — Mquj‘ t sincce 


dz ] -(- c.c. , 


( 5 ) 


in the Cartesian coordinates associated with the Boyer- 
Lindquist frame, with the upper (lower) sign correspond- 









3 


ing to an aligned (anti-aligned) binary system. We will 
also consider the case where any deformation rotates with 
the pulsar, such that its quadrupole moment corresponds 
to that of a deformed ellipsoid with principal axis of in¬ 
ertia at an angle (3 from the spin axis. The resulting 
quadrupole tensor includes two frequency components, 
as anticipated above, yielding [11]: 


Qij 

Qo 



0 0 ±i 
0 0-1 
±i -1 0 


+ c.c. 


/I ±i 0\ 

+ spe~ 2iuipt l ±i -1 0 

\ 0 0 0 / 

( 6 ) 


where cp = cos/3 and sp = sin j3. The amplitudes Mq and 
Q o can be obtained from the properties of the neutron 
star but will not be relevant to our discussion. We note 
that the complex conjugate terms in Eqs. ( 0 ) and © lead 
to waves of the form e lui . 

From this we can use the standard dipole and 
quadrupole formulas to determine the EM and GW 
fields generated by the pulsar and the corresponding NP 
scalars. In the pulsar’s frame, the s < 0 NP scalars are, 
up to an s— dependent magnitude, given by: 

-iu(t-r) „iuj(t-r) 

T s --- s Y lp:mp (0 » + ^— s Y lp ,-m P {0A) ,(7) 


where l p = 1 and m p = ±1 for EM radiation of frequency 
uj p , while for GWs we have l p = 2 and m p = ±1,±2 for 
u) = Lu p , 2w p , respectively. 

In the BH frame, however, the NP scalars become 
plane waves for large orbital distances, which are super¬ 
positions of all the different multipoles l > |s|. The mul¬ 
tipolar decomposition of the NP scalars, which we will 
describe in detail in a companion article [171 ]. yields a 
remarkably simple result for the incoming radiation (di¬ 
rected towards the BH). For the e~ lLJt modes, we obtain 
for the incident power: 


dE). 


(s,Z,m) 


dt 


= ^~d 


sY lp ^ mp {0 p )\ 2 \ s Y h _ m (9 p )\ 2 P s (8) 


where P s denotes the pulsar’s total luminosity in EM 
or GWs, while for the e lLJt modes we find an analogous 
result with m p —> —m p . 

We thus conclude that the fraction of the incident 
power carried by each (l,m) mode depends on the po¬ 
lar angle 9 p but not on the azimuthal angle <f > p , which 
only determines the phase of each mode. To better illus¬ 
trate this, we show in Fig. [L] the relative contribution of 
the different GW quadrupole modes. 

As shown in this figure, the m > 0 modes are dominant 
at large angles, while the m < 0 modes give the largest 
contribution at small angles. Close to the BH’s equator, 
9 P = 7t/2, most of the radiation is in the non-superradiant 
m = 0 mode, and modes of opposite m are present in 
equal amounts. This is a generic behavior, common to 
all multipoles in both the GW and EM cases. 


0 2 7T 
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FIG. 1: Relative fraction of the incident energy carried by 
the l = 2 modes of the gravitational radiation emitted by the 
pulsar, normalized to the total quadrupole contribution, as a 
function of its polar angular position, 6 P . 


Additionally, the relative contribution of the e ±KlJt 
modes also depends on the angle 9 p through distinct spin- 
s harmonics, which correspond to the NP scalars in the 
pulsar’s frame. Hence, the amount of (non-)superradiant 
modes in the incident radiation depends on the pulsar’s 
angular position relative to the BH’s rotation axis, and 
the radiation is generically “polarized”. 

The fraction of the incident energy that is absorbed or 
amplified by the BH can be determined by multiplying 
the incident energy of each mode by the corresponding 
gain/loss factor s Zi m . We may then obtain the effective 
pulsar luminosity by subtracting the energy flowing into 
the BH per unit time from the total power. We illustrate 
this in Fig. [2] for GWs of frequency uj = 2 uj p ~ 2f 1 and 
a = 0.999. This would correspond to e.g. a millisecond 
pulsar around a 16 Mq near-extremal BH. 

This figure shows that results depend on the relative 
orientation of the pulsar and BH spins. On the one hand, 
when the spins are aligned, superradiant modes of the 
form e~ lult and m > 0 ( e lut and m < 0) are domi¬ 
nant at large (small) angles, and the signal is amplified 
when the pulsar’s position makes an angle larger than 
~ 62° with the equatorial plane, both above and below 
the equator. Closer to the equator, the incident radiation 
is dominated by non-superradiant modes and the pulsar 
luminosity decreases. We note that, in this example, only 
the l = 2,3 modes are relevant, with higher multipoles 
having negligible gain/loss factors. 

On the other hand, when the spins are anti-aligned, 
the relative contributions of the and e lult modes are 
reversed. The incident wave is then dominated by non¬ 
superradiant modes and the radiation is mostly absorbed 
for all angles. In the equatorial plane, modes with oppo¬ 
site m and w are present in equal amounts and the signal 
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FIG. 2: Modulation of the pulsar’s luminosity in GWs of fre¬ 
quency w = 2 cjp ~ 2fi with a near-extremal BH companion, 
a = 0.999, as a function of the polar angle 9 P . The solid blue 
(dashed green) curve corresponds to the case where the pulsar 
and BH spins are aligned (anti-aligned). The orbital distance 
is d = 10A in this example. 



FIG. 3: Modulation of a pulsar’s GW luminosity for a polar 
circular orbit around a near-extremal BH (5 = 0.999), for fre¬ 
quency ui = 2 cj p ~ 2ft and orbital radius d = 10A. The solid 
blue (dashed green) curve corresponds to the case where the 
pulsar and BH spins are aligned (anti-aligned). The bottom 
diagrams depict the pulsar’s position in its orbit, with the 
arrow indicating the BH spin axis. 


is attenuated by the same factor in both the aligned and 
anti-aligned configurations. 

We find similar results for EM waves, although amplifi¬ 
cation is less pronounced (< 10 -4 ) and requires angular 
deviations from the equatorial plane of more than 77°. 
The example in Fig. [2] corresponds to the most promising 
scenario, since amplification factors are maximal for GWs 
close to the superradiance threshold ui ~ 2fl and near¬ 
extremal BHs. For lower frequencies and BH spins we 
nevertheless find similar qualitative results, with e.g. an 
overall modulation of ~ 0.01% for a = 0.9, for which we 
will give a detailed analysis in 171. 

Orbits in the Kerr space-time are generically non- 
planar, making the pulsar probe different polar angles 
in its orbit [19j and thus inducing a modulation of its 
luminosity. Here we consider the simplest case of pla¬ 
nar circular orbits, which are a good approximation at 
large orbital radius. Since the effects of superradiance 
are more pronounced close to the poles, we illustrate in 
Fig. [3] the best-case scenario of a polar orbit for the ex¬ 
ample considered above. 

While we observe a single modulation of the luminosity 
for anti-aligned spins, an aligned binary system exhibits 
an additional modulation, which is entirely due to super¬ 
radiant scattering when the pulsar moves away from the 
equatorial plane. This double-modulation is a generic fea¬ 
ture, which is more or less pronounced depending on the 
radiation frequency, BH spin and orbital inclination, thus 
potentially yielding a ‘smoking gun’ for superradiance. 

Although we have considered the total pulsar luminos¬ 
ity, a similar modulation will occur along different lines- 
of-sight, as we will describe in detail in 117]. Note that in¬ 
coming modes are further modulated by the correspond¬ 
ing spin-s harmonics along a given direction, which may 


either enhance or reduce the signal modulation at differ¬ 
ent points in the orbit. The analysis given above never¬ 
theless yields the generic magnitude and orbital evolution 
of the effects of superradiance. Weak lensing may also po¬ 
tentially pose an important experimental challenge along 
particular lines-of-sight 0], although its effects may in 
principle be distinguished from superradiant amplifica¬ 
tion by measuring different EM and GW frequencies. 

In the example above the effects of superradiance on 
the pulsar’s luminosity are below the percent level. How¬ 
ever, we expect this effect to become more pronounced for 
smaller orbital radii, namely for systems close to merg¬ 
ing, where the incident power is larger. Moreover, the 
incident wave should approach the single multipole form 
obtained in the pulsar’s frame as the radius decreases, 
thus favoring superradiant modes for an aligned binary. 
Although this lies outside the scope of this work, it mo¬ 
tivates further study of this problem beyond the plane 
wave approximation. 

For EM waves, the modulation effect is quite sup¬ 
pressed and, furthermore, low-frequency radio waves are 
easily absorbed/attenuated in astrophysical plasmas. In 
fact, the plasma in the pulsar’s magnetosphere may pre¬ 
vent the magnetic dipole radiation from reaching the BH. 
Even for pulsars where this plasma is periodically absent, 
such as B1931+24 [HI, free-free absorption in the inter¬ 
stellar medium will most likely prevent the radiation from 
reaching the Solar System, and the superradiant modu¬ 
lation can at most be probed through its indirect effect 
e.g. in heating intervening gas clouds 22|. 

GWs can yield a much stronger and cleaner evidence 
for superradiance in pulsar-BH binaries. The next gen¬ 
eration of ground-based detectors, such as Advanced 
LIGO/VIRGO, should be sensitive to GWs emitted by 
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deformed neutron stars, with current detectors alr eady 
placing stringent constraints on pulsar ellipticities [231 ]. 
Our results thus give additional motivation to search for 
pulsar gravitational radiation in such binary systems. 

Evidence for BH superradiance would be a very impor¬ 
tant step in testing the yet unexplored corners of general 
relativity. Superradiance may also yield a unique labora¬ 
tory for particle physics, through the potential instabili¬ 
ties that it may generate for ultra-light fields confined in 
the BH’s gravitational potential ]3, 2J-36|. This can also 
lead to massive scalar ‘hair’ around a Kerr BH [37]. 

The pulsar luminosity modulation found in this Let¬ 
ter provides the first of potentially other ways of testing 
the occurrence of superradiant scattering in astrophysical 
BHs and we hope that this motivates further exploration 
of this fascinating aspect of BH physics. 
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